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ABSTRACT 
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conducted.  The  recursive  ray  acoustics  (RRA)  algorithm  was  found  to  be  accurate  and 
relatively  fast.  The  RRA  algorithm  is  capable  of  handling  sound  speed  as  a  function  of 
all  three  spatial  coordinates,  and  this  capability  was  demonstrated. 

Two  separate  methods  of  representing  a  sound-speed  profile  (SSP)  based  on  data 
points  were  examined:  Akima  cubic  spline  and  spatial  Fourier  series  (SFS).  The  SFS 
representation  encountered  difficulties  in  accurately  modeling  SSPs.  Various  techniques 
were  applied  to  improve  the  SFS  sound-speed  representation.  While  accurate  sound- 
speed  fits  were  eventually  achieved,  difficulties  remained  in  the  SFS  modeling  of  first  and 
second-order  derivatives  of  the  sound-speed  data. 

The  RRA  algorithm  was  tested  using  the  SFS  sound-speed  representation  and  found 
to  be  significantly  inaccurate. 

A  demonstration  was  conducted  of  the  ability  of  the  SFS  sound-speed  represen- 
tation to  incorporate  randomness  in  the  SSP. 
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I.     INTRODUCTION 

The  primary  focus  of  this  thesis  was  to  test  the  accuracy  and  capabilities  of  a  re- 
cursive ray  acoustics  (RRA)  algorithm  that  can  be  used  to  calculate  the  position,  angles 
of  propagation,  travel  time,  and  path  length  along  a  ray  path  and  to  generate  ray  trace 
plots  for  speeds  of  sound  that  are  functions  of  all  three  spatial  coordinates.  The  RRA 
algorithm  is  an  extension  and  slight  modification  of  an  algorithm  briefly  discussed  by 
Klein  [Ref.  1:  pp.  29-31].  The  algorithm  discussed  by  Klein  [Ref.  1:  pp.  29-31  ]  was  ex- 
pressed in  terms  of  an  index  of  refraction  while  herein,  the  RRA  algorithm  is  expressed 
directly  in  terms  of  the  speed  of  sound.  No  computer  simulation  results  were  presented 
by  Klein  [Ref.  1:  pp.  29-31  ]. 

The  accuracy  of  the  RRA  algorithm  was  tested  by  comparing  its  computer  simu- 
lation results  with  those  obtained  from  a  previously  developed  ray  acoustics  algorithm 
based  on  solving  a  system  of  first-order  ordinary  differential  equations  (ODE)  [Ref.  2]. 
The  ODE  algorithm  as  applied  to  ray  acoustics  was  discussed  in  detail  by  Lim  [Ref.  3]. 
Both  the  ODE  algorithm  and  the  RRA  algorithm  used  the  IMSL  Version  10  double 
precision  computer  program  DCS  A  KM  to  obtain  an  Akima  cubic  spline  fit  to  one- 
dimensional,  depth-dependent,  sound-speed  data  for  testing  and  comparison  purposes. 
The  computer  simulation  results  are  presented  in  Chapter  III,  Section  A.  An  Akima 
cubic  spline  was  chosen  since  it  is  designed  to  produce  a  curve  that  matches  the  shape 
of  the  data  while  minimizing  oscillations.  In  addition,  the  RRA  algorithm  was  tested 
and  shown  to  be  capable  of  handling  a  sound-speed  profile  (SSP)  that  is  a  function  of 
all  three  spatial  coordinates.  This  capability  was  demonstrated  by  first  separating  the 
ocean  into  zones  along  the  z  axis.  In  each  zone,  the  SSP  is  a  different  function,  in  gen- 
eral, of  both  cross-range  x  and  depth  y.  A  SSP  with  a  gradient  in  the  cross-range  di- 
rection will  cause  a  sound  ray  to  leave  its  initial  plane  of  propagation.  Chapter  III, 
Section  B,  includes  demonstrations  of  the  special  three-dimensional  capabilities  of  the 
RRA  algorithm. 

A  secondary  objective  of  this  thesis  was  to  examine  the  use  of  a  spatial  Fourier  series 
(SFS)  to  fit  sound-speed  data  and  to  assess  its  utility  relative  to  an  Akima  cubic  spline 
fit.  The  appeal  of  the  SFS  lay  in  its  representation  of  spectral  content  in  terms  of  spatial 
frequencies,  as  well  as  its  ability  to  represent  derivatives.  The  precision  of  the  two 
sound-speed  fits,  as  well  as  that  of  the  first  and  second-order  derivatives  of  the  sound- 


speed  data  were  compared  since  the  first-order  derivative  as  well  as  the  sound-speed  fit 
itself  are  utilized  in  the  ray  trace  algorithms  mentioned  above.  Chapter  IV,  Section  A, 
contains  the  results  of  this  comparison  as  generated  by  computer  simulation.  The  SFS 
representation  of  a  SSP  was  then  employed  within  the  RRA  algorithm  and  tested  by 
using  the  same  SSPs  that  were  used  to  compare  the  ODE  and  RRA  algorithms.  Chapter 
IV,  Section  B,  discusses  the  accuracy  of  the  RRA  algorithm  incorporating  the  SFS  rep- 
resentation of  the  various  SSPs. 

Additionally,  random  fluctuations  in  the  SSP  were  modeled  using  the  inherent  suit- 
ability of  a  SFS  representation  for  the  inclusion  of  components  at  particular  harmonics. 
This  random  element  was  intended  to  simulate  fluctuations  generated  by  such  processes 
as  internal  waves  or  by  inaccurate  collection  of  sound-speed  data. 


II.     THEORETICAL  ANALYSIS 

A.     THE  RRA  ALGORITHM 

If  a  sound  source  is  located  at  r0  =  {x0,y0,  zo)  and  if  a  ray  is  launched  with  initial  an- 
gles of  propagation  P0  and  0O  (see  Figure  1  ),  then  an  estimate  of  the  position  of  subse- 
quent points  along  the  ray  path  is  given  by  (see  Figure  2  ) 

4  =  4-1  +&Zi-v    /=  1,2,3, ...  (2.1) 

where 

A^  =  A5i«/,  (2.2) 

As,  is  an  increment  of  arc  length,  the  magnitude  of  which  is  at  the  user's  discretion, 

nt  =  U[X  +  vty  +  wtz  (2.3) 

is  the  unit  vector  along  the  ray  path  where 

u,  =  sin  Pi  cos  4>i ,  (2.4) 

v^cos  p^  (2.5) 

and 

wt  =  sin  Pi  sin  #,  (2.6) 

are  the  dimensionless  direction  cosines  with  respect  to  the  x,y  and  z  axes,  respectively, 
and  P,  and  </>,  are  the  angles  of  propagation. 

Since,  in  general,  «,_,  is  not  identical  with  h„  from  Klein  [Ref.  1:  pp.  29-31]  we  see 
that,  as  amplified  and  applied  to  the  ocean  acoustics  environment  in  Ziomek  [Ref.  4:  pp. 
225  and  230], 

d(n(r)n{r)) 

ds         ~WnU)'  (2,?) 


Figure   1.      Initial  Angles  of  Propagation  /?0  and  <f>0 
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Figure  2.      Ray  Geometry. 


or,  approximately, 

A(nfo$)  -  Vn{^ASl  (2.8) 

where  h(r)  =  h, , 

c(fy) 

"W=W  <2-9> 

is  the  dimensionless  index  of  refraction,  c(£,)  is  the  three-dimensional  speed  of  sound  at 

rQ  =  xQx+yQy  +zQz  ,  (2.10) 

and,  c(xj)  is  the  three-dimensional  speed  of  sound  at 

Zi  =  xlx+y1y  +ztz  ,  (2.11) 

which  is  the  position  vector  to  a  point  on  the  ray  path  after  the  ith  iteration. 
It  follows  from  Equation  (2.8)  that 

nfe)  nt  -  nir^n^  =  Vn{ri)Asi  (2.12) 

and,  as  a  result, 

A      «(i:,_t  )«,•_,       As i         * 

where,  for  better  accuracy  in  the  evaluation  of  the  gradient, 

*  ASl-\       A 

r,_!  =  r,_,  +  — 2~  «/_i.  (2-14) 

which  is  the  three-dimensional  spatial  location  half-way  along  the  arc  length  between 
£t_i  and  r.i,  and 

V  =-r—  *  +-^—  y  +-r~  z  (2-15) 

3jc  dy J      dz  v 

is  the  gradient  expressed  in  rectangular  coordinates.    It  should  be  noted  that  although 


As,  has  been  shown  as  a  function  of  iteration  number,  it  will  be  constant,  in  general, 
except  when  it  needs  to  be  adjusted  to  properly  handle  surface  and  bottom  reflections 
and  maximum  down  range  limits. 

In  terms  of  the  direction  cosines   u,  v,  and  w, 


«(£_ ,)«,_,         As     Mti-i) 


n(n)  /•?(/,)        dx 

«Gi-ih-i         As     d«G£-i) 


+ 


n{Lt)  «fo)        dy 


(2.16) 
(2.17) 


and 


+ 


As     Md-i) 
nfe)        dz 


(2.18) 


If  c(r,)  =  c(jc„j;„  zj  is  the  three-dimensional  speed  of  sound,  and  if  c0  =  c^)  is  the 
speed  of  sound  at  the  source  location,  then  it  follows  that 


"(r(_,)  c0       cfo)         cfo) 


«(ii)         cfo.j)     co         cfo^) 


(2.19) 


and 


1       dnir) 
n{r-)      dx 


=  ~  cfo) 


l      3c(r) 


c(L) 


dx 


(2.20) 


Analogous  reasoning  applied  to  the  y  and  z  coordinates  yields  results  similar  to 
Equation  (2.20)  ,  and  substitution  of  these  results  into  Equation  (2.13)  yields 


cfo-i) 


«/_]  -  cfo)As, 


1 

L  c2(.) 


Vc(r) 


'/-i 


(2.21) 


In  addition,  although  Klein  did  not  discuss  travel  time  t  and  path  length  5  calcu- 
lations, these  quantities  were  computed  using  the  following  formulae: 


r=    >  —+-,  (2.22) 


and 


-I 


As, .  (2.23) 

Equations  (2.1)  through  (2.23)  constitute  what  is  referred  to  herein  as  the  recursive 
ray  acoustics  (RRA)  algorithm  first  mentioned  in  Chapter  1.  Input  to  this  algorithm  in 
the  form  of  mathematical  formulae  for  the  sound-speed  profile  or  some  particular  nu- 
merical fit  to  sound-speed  data  must  be  provided  to  enable  the  algorithm  to  function. 

B.     SPATIAL  FOURIER  SERIES  REPRESENTATION  OF  SOUND-SPEED  DATA 

Given  an  ocean  D  meters  deep  with  a  one-dimensional,  depth-dependent  sound- 
speed  profile  (SSP),  the  speed  of  sound  c(y)  can  be  represented  by  the  following  spatial 
Fourier  series  (SFS): 

5 

c(y)  =  I  c0  I  +  2^  I  cs\  cos(2nfyoy  -  Lcs),  0<y<D  (2.24) 

5=1 

where  |  cs  |  and  Lcs  are  the  magnitude  and  phase(angle),  respectively,  of  the  complex 
Fourier  series  coefficient  c,  at  harmonic  s,  \c0\  is  the  magnitude  of  cs  at  harmonic  5  =  0 
and  is  the  "DC  component"  or  average  value  of  the  sound-speed  data, 

/  =  — L-  (2.25) 

is  the  fundamental  spatial  frequency  in  the  j>  direction  in  cycles  per  meter  where  N  is  the 
total  number  of  sound-speed  data  points  taken  and  Ts  is  the  uniform  sampling  period 
in  the  y  direction  in  meters  per  sample,  and  S  is  the  highest  harmonic.  Derivatives  of 
this  Fourier  series  representation  can  be  directly  obtained. 


The  complex  Fourier  series  coefficients  cs  are  obtained  by  taking  the  spatial  discrete 
Fourier  transform  (DFT)  of  the  sound-speed  data  as  follows: 

<v-i 
cs—fi-Yie(nTs)W"'    *--S,.»,Of...fS  (2.26) 

where 

lV„=exp(+j7nlN)  (2.27) 

and,  when  sampling  at  the  Nyquist  rate, 

N=2S+l.  (2.28) 

An  attempt  to  represent  a  function  by  a  truncated  Fourier  series  will  always  result 
in  oscillations  about  the  functions'  true  value.  These  oscillations  normally  have  small 
enough  amplitudes  that  they  are  of  little  consequence.  However,  an  attempt  to  use  a 
SFS  representation  of  a  linear  SSP  with  a  single  constant  gradient  results  in  significant 
discrepancies  in  amplitude  near  the  ocean  surface  and  bottom,  due  to  the  inherent  dis- 
continuity in  the  SSP  that  exists  at  these  boundaries.  This  difficulty,  known  as  "Gibbs' 
phenomenon",  does  not  diminish  as  more  harmonics,  that  is,  sound-speed  data  points, 
are  added  to  the  Fourier  series. 

C.  Lanczos  [Ref.  5]  suggested  that  this  phenomenon  would  be  ameliorated  if  the 
Fourier  series  could  be  forced  to  converge  more  rapidly  by  modifying  the  Fourier  series 
coefficients  cr  by  a  correction  factor,  such  that  the  Lanczos  smoothed  Fourier  series  co- 
efficients cSL  are  given  by 

sin(s7r/S) 

^-Hbs    *  *#0'  (229) 

where  c,  is  given  by  Equation  (2.25).  Hamming  [Ref.  6]  states  that  S  in  Equation  (2.29) 
is  usually  replaced  with  S  -1-  1  and,  thus,  our  representation  of  the  correction  factor  will 
be 

sin(57t/(S+l))  ,„^M 


C.     INCORPORATION  OF  RANDOMNESS  INTO  THE  SOUND-SPEED 
PROFILE 

The  speed  of  sound  as  a  function  of  depth,  c(y)  ,  can  be  represented  as  the  sum  of 
deterministic  and  random  components  as  follows: 

c(y)  =  cD(y)  +  cR(y),  (2.31) 

with  cD(y)  being  the  deterministic  and  cR(y)  the  random  speeds  of  sound  as  functions  of 
depth y. 

The  autocorrelation  function  of  cR(y)  is  ,  for  any  particular  y  andy  , 

RCRiy,y')  =  E{cR{y)c*R{y')),  (2.32) 

where  cR{y')  is  the  complex  conjugate  of  cR  at  depth  y,  0  <,y'  <,  D  . 
If 

N 

ci&)  =  X  Cn  exp(  ~  J2™fy*y)>  (2-33) 

for  some  set  of  Fourier  series  coefficients  cn,  and 

4=1/0  (2-34) 

for  D  the  ocean  depth  in  meters,  then 

1    tD 

Cn  =  ~o)  crty)QXP(+J2nnfy0yWy-  (2-35) 

If  cR(y)  is  a  random  function  of  depth,  then  cn  will  be  a  random  variable. 
Since 


%»(<V)*}  =  A"  \D{DE{cR(y)cR(y')}  exp(  +  jlnfyo(ny  -  n'y'))dydy'  ,  (2.36) 


where 

E{cR(y)cR(y')}  =  RCR{y,y'),  (2.37) 


and  (c„,)'  is  the  complex  conjugate  of  cn. ,  and  if  it  is  assumed  that  cR{y)  is  a  wide-sense 
stationary  (  i.e.,  statistically  homogeneous  )  random  process,  then 

RCR(y,y')  =  RCR(Ay),  (2.38) 

where  Ay  =  y  —  /  .  If  R     is  represented  by  a  finite  Fourier  series 

N 
RCR(Ay)  =   Y,  rm  exp(  -  J2nmfyAy) ,  (2.39) 

then,  as  developed  in  Thomas  [Ref.  7:  pp.  148-153], 

r™  =  TnT\  s^mfy^  (2-40) 

where  SCR  is  the  power  spectral  density  of  cR{y). 

Thomas  [Ref.  7:  pp  148-153]  goes  on  to  show  that 

E{cnc»]  =  rnSnn,,  (2.41) 

where  Sm,  is  the  Kronecker  delta  function.   Therefore, 

Hence,  if  cn  is  a  zero  mean  Gaussian  random  variable  with  variance  rn  ,  that  is,  if 

cn  =  N(0,  rn),  (2.43) 

then  cR(y)  will  be  a  zero  mean  Gaussian  random  process. 
Therefore, 

Rcjy.y)  =  *4*v)  =  Er"  exp(  -  J2nfyM  (2-44) 

n 

which  represents  the  Inverse  Discrete  Fourier  Transform  of  the  power  spectral  density. 
Thus,  knowledge  of  the  power  spectral  density  of  cR(y)  enables  calculation  of  rn  via 
Equation  (2.42),  which  in  turn  yields  cn  via  (2.43)  and  cR(y)  via  (2.33). 
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III.     COMPUTER  SIMULATION  RESULTS 

A.     RAY  ACOUSTICS  -  COMPARATIVE  ANALYSIS  OF  ODE  VERSUS  RRA 
ALGORITHMS 

Four  distinct  sound-speed  profiles  (SSPs)  were  examined  for  purposes  of  compar- 
ative analysis.  For  all  cases,  a  sound  source  spatial  location  of  x0  =  0  m  and  z0  =  0  m 
was  used.  The  source  depth  y0  was  established  at  a  point  half-way  between  the  ocean 
surface  and  bottom,  100  m  in  most  cases.  A  constant  ocean  depth  of  200  m  was  mod- 
eled, unless  otherwise  noted,  and  for  all  cases  three  separate  rays  were  examined.  All 
rays  were  propagated  initially  in  the  direction  of  the  positive  z  axis,  that  is,  0O  =  90°,  but 
/?0  was  altered  in  order  to  generate  unique  rays.  For  SSPs  1  to  3,  the  launch  angle  /?0 
was  set  to  45°,  85°,  and  135°  .  SSP  4  was  designed  to  simulate  a  SOFAR  channel. 
Launch  angles  were  therefore  chosen  to  result  in  the  first  two  rays  propagating  down 
range  without  reflection,  i.e.,  /?0  =  85°  and  95°,  while  the  third  launch  angle  was  set  to 
/?0  =  135°  to  test  reflection  as  well. 

The  four  SSPs  considered  are  defined  as  follows: 

•  SSP  1:  Constant  speed  of  sound  (  c  =  1500  m/s  ) 

•  SSP  2:  Linear  SSP  with  a  positive  gradient  (  c(y)  =  c(0)  +  gy  m/s  where  the  speed 
of  sound  at  the  surface  c(0)  =  ciy  =  0)  =  1500  m/s  and  g  =    +0.017  1/s  ) 

•  SSP  3:  Linear  SSP  with  a  negative  gradient  (  c{y)  =  c(0)  +  gy    m/s  where  the 
speed  of  sound  at  the  surface  c(0)  =  c(y  =  0)  =  1500  m/s  and  g  =    —0.017   1/s  ) 

•  SSP  4:  Parabolic  SSP  with  vertex  at  y0  where  c(y0)  =  1490  m/s  and  the  speed  of 
sound    at    both    the    ocean    surface    and    bottom    is    equal    to     1500    m/s 

(  c(y)  =  c{y0)  +  a{y  —  y0)2    m/s  with  a  =  0.001    -^^r  for  an  ocean  depth  of  200 
m). 

The  above  four  SSPs  generated  nine  sound-speed  data  points  which  were  input  to 
the  Akima  cubic  spline  routine.  The  results  from  the  ODE  algorithm  were  used 
throughout  as  the  benchmark  for  comparison  purposes.  The  ODE  algorithm  used  the 
IMSL  Version  10  double  precision  computer  program  DIVPBS  in  order  to  solve  the 
system  of  first-order  ordinary  differential  equations.  A  range  step  size  of  2.0  m  and  an 
error  tolerance  of  1  x  10~7  was  used  in  DIVPBS,  except  for  SSP  3,  Ray  1  in  Table  3, 
where  the  tolerance  had  to  be  reduced  to  1  x  10-6  in  order  for  the  ODE  algorithm  to  run. 
The  arc  length  step  size  used  in  the  RRA  algorithm  was  2.0  m  unless  otherwise  noted. 
The  results  of  the  various  simulation  runs  are  presented  in  the  following  seven  tables. 
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Table  1.     RAY  ACOUSTICS  RESULTS  FOR  SSP  1,  RAYS  1  THROUGH  3 


Method 

}'o  (m) 

&  (dcg) 

r  (km) 

y  (m) 

P  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

100.00 

45.000 

9.428090 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.95 

45.000 

9.428091 

14.1421 

ODE 

100.0 

85.0 

10.0 

174.89 

85.000 

6.692132 

10.0382 

RRA 

100.0 

85.0 

10.0 

174.88 

85.000 

6.692132 

10.0382 

ODE 

100.0 

135.0 

10.0 

100.00 

135.000 

9.428090 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.05 

135.000 

9.428091 

14.1421 

For  SSP  1  (constant  speed  of  sound),  Table  1  shows  that  at  a  horizontal  range 
r  =  10  km,  the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.01  and  0.05  m  in  depths  and 

•  1  ix  sec  in  travel  time  r. 
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Table  2.     RAY  ACOUSTICS  RESULTS  FOR  SSP  2,  RAYS  1  THROUGH  3 


Method 

y*  (m) 

Po  (deg) 

r  (km) 

y(m) 

fi  (deg) 

t  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.417413 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.95 

45.000 

9.417422 

14.1421 

ODE 

100.0 

85.0 

10.0 

162.50 

85.487 

6.683492 

10.0374 

RRA 

100.0 

85.0 

10.0 

162.57 

85.487 

6.683496 

10.0374 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.417413 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.05 

135.000 

9.417422 

14.1421 

For  SSP  2  (linear  SSP  with  a  positive  gradient),  Table  2  shows  that  at  a  horizontal 
range  r  =  10  km,  the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.04  and  0.07  m  in  depth  y  and 

•  between  4  and  9  (x  sec  in  travel  time  x. 
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Table  3.     RAY  ACOUSTICS  RESULTS  FOR  SSP  3,  RAYS  1  THROUGH  3 


Method 

Vo  (m) 

A  (deg) 

r  (km) 

y  (m) 

fi  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.438784 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.95 

45.000 

9.438795 

14.1421 

ODE 

100.0 

85.0 

10.0 

168.89 

84.512 

6.699404 

10.0380 

RRA 

100.0 

85.0 

10.0 

168.96 

84.511 

6.699409 

10.0380 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.438784 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.05 

135.000 

9.438794 

14.1421 

For  SSP  3  (linear  SSP  with  a  negative  gradient),  Table  3  shows  that  at  a  horizontal 
range  r  =  10  km,  the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.04  and  0.07  m  in  depth  y  , 

•  0.001  deg  in  angle  of  arrival  /?  ,  and 

•  between  5  and  1 1  \x  sec  in  travel  time  r. 

The  data  presented  in  Tables  1  through  3  demonstrate  that  the  two  algorithms  are 
in  very  good  agreement.  Note  that  the  accuracy  of  the  RRA  algorithm  can  always  be 
improved  by  reducing  the  arc  length  step  size.   This  is  demonstrated  in  Tables  5  and  6. 
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Table  4.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3 


Method 

>'o  (m) 

&  (deg) 

r  (km) 

y  (m) 

P  (deg) 

t  (sec) 

s  (km) 

ODE 

100.0 

85.0 

10.0 

37.06 

87.244 

6.710387 

10.0183 

RRA 

100.0 

85.0 

10.0 

36.63 

87.224 

6.710335 

10.0184 

ODE 

100.0 

95.0 

10.0 

162.94 

92.756 

6.710387 

10.0183 

RRA 

100.0 

95.0 

10.0 

163.37 

92.776 

6.710335 

10.0184 

ODE 

100.0 

135.0 

10.0 

144.9^ 

134.922 

9.448852 

14.1104 

RRA 

100.0 

135.0 

10.0 

144.24 

134.930 

9.449238 

14.1110 

For  SSP  4  (parabolic  SSP),  Table  4  shows  that  at  a  horizontal  range  r  =  10  km,  with 
the  ocean  bottom  at  200  m,  the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.43  and  0.73  m  in  depth  y  , 

•  between  0.008  and  0.02  deg  in  angle  of  arrival  /?  , 

•  between  52  and  386  ju  sec  in  travel  time  t  ,  and 

•  between  0.1  and  0.6  m  in  path  length  s. 

As  is  indicated  in  Table  4,  the  differences  in  travel  time  are  relatively  large  when 
compared  with  previous  cases.  This  lower  level  of  accuracy  was  considered  unacceptable 
since,  in  general,  accuracy  in  the  millisecond  region  is  normally  required  for  travel  time. 
An  attempt  to  improve  accuracy  was  made  by  reducing  the  arc  length  increment  used 
in  the  RRA  algorithm  from  2.0  to  1.0  m  (see  Tables  5  and  6). 

The  parabolic  SSP  which  produced  Table  4  has  the  most  extreme  gradient  of  all  the 
SSPs  examined  and,  thus,  the  two  rays  launched  at  85°  and  95°  proceed  down  range 
without  reflection  from  the  surface  or  bottom,  as  illustrated  in  Figures  3  and  4.  This 
extreme  gradient  was  considered  to  have  been  the  probable  cause  of  the  discrepancies 
noted  above  and,  therefore,  the  gradient  was  reduced  and  the  ocean  depth  increased  to 
permit  an  examination  of  a  similar  profile  with  a  more  typical  gradient  (see  Table  7). 


15 


o 

in    i 

t^~f 

CO 

o 

1 — 1 

en 

ID 
O 
CD 

en 

>-« 

az 
en 

> 
• — i 

CO 

en 

ID 
CJ 
L±J 
C£ 

in 

1 

? 

X 

1      "•_/ 
i     **    I 

in  i 

CD 

UJ 

o 
o 

o  .. 

™2 

cr> 

-»-  t"1  CD 

a:  2^0 

ceoLO 
•  to 
v  o 
o  .. 

cj      az 

CO   ••  1— 
CC  OLJ 
CJ>->0Q 

s 

8 

C3 

ft 

O    1 

Figure  3.      SSP  4  -  No  Reflections  for  /J0  =  85 


16 


o 

„ 

in    I 

cs>"  ~f 

i*:  ir>- 

■ 

C*J 

CO 

o 

1 — 1 

in        / 

I      1                1 

C9 

LJ 

o 

CO 
ZD 
O 

iW 

o 

_  D 
31  (J) 

O 

cn 

O  .. 

•o 

CM.-! 

>- 

az 
cn 

LJ 

> 

t— 1 
CO 

a: 

Zj 
CJ 
LJ 
C£ 

*■*    y 

3: 
l>Jcn 

CC3=0 

ctroin 
•  cn 

o  .. 
••  —  a 
u      cc 

CO   ■■  H> 

ccou 

0>-tCQ 

2_;      J 

X         \ 

3  1 

8 

8 

a 

% 

o    1 

Figure  4.      SSP  4  -  No  Reflections  for  /J0  =  95 


17 


Table  5.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

Jo  (m) 

Po  (deg) 

r  (km) 

y(m) 

fi  (deg) 

t  (sec) 

s  (km) 

ODE 

100.0 

85.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

64.13 

85.581 

3.354757 

5.0088 

ODE 

100.0 

95.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

95.0 

5.0 

135.87 

94.419 

3.354757 

5.0088 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

135.0 

5.0 

77.66 

45.017 

4.724506 

7.0553 

The  data  in  Table  5  for  a  horizontal  range  of  5  km  and  using  an  arc  length  step  size 
of  2.0  m  was  generated  in  order  to  enable  a  direct  comparison  with  the  arc  length  step 
size  in  the  RRA  algorithm  set  to  1.0  m,  as  shown  in  Table  6.  The  reduction  in  final 
range  from  10  to  5  km  was  necessitated  by  the  structure  of  the  RRA  computer  code, 
which  required  large  quantities  of  computer  storage  space  to  be  allocated,  particularly 
as  arc  length  step  size  was  decreased.  This  drawback  could  be  rectified  by  straightfor- 
ward modifications  to  the  simulation  program. 

For  SSP  4  (parabolic  SSP),  Table  5  shows  that  at  a  horizontal  range  of  r  =  5  km, 
the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.12  and  0.16  m  in  depth  y  , 

•  between  0.002  and  0.016  deg  in  angle  of  arrival  /?  , 

•  between  22  and  95  n  sec  in  travel  time  x  ,  and 

•  0.1  m  in  path  length  s 
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Table  6.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  1.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

Vo  (m) 

0o (deg) 

/•  (km) 

y(m), 

fi  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

85.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

64.19 

85.589 

3.354768 

5.0088 

ODE 

100.0 

95.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

95.0 

5.0 

135.81 

94.411 

3.354768 

5.0088 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

135.0 

5.0 

77.56 

45.018 

4.724458 

7.0553 

For  SSP  4  (parabolic  SSP),  Table  6  shows  that  at  a  horizontal  range  of  r  =  5  km, 
and  with  the  RRA  arc  length  step  size  set  to  1.0  m,  the  differences  between  the  two  al- 
gorithms are  as  follows: 

•  0.06  m  in  depth  y  , 

•  between  0.001  and  0.008  deg  in  angle  of  arrival  /?  , 

•  between  1 1  and  47  /j.  sec  in  travel  time  t  ,  and 

•  0.1  m  in  path  length  s. 

A  significant  improvement  in  the  consistency  of  travel  time  calculations  can  be  seen 
by  comparing  Tables  5  and  6.  Cutting  the  arc  length  increment  in  half  reduced  the  travel 
time  discrepancy  by  approximately  one-half.  Further  reductions  in  arc  length  increment, 
an  input  parameter  to  the  simulation  program,  would  yield  further  improvements  in 
travel  time  accuracy. 
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Table  7.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  BOTTOM  DEPTH  1000  M 


Method 

y0  (m) 

P0  (deg) 

r  (km) 

y(m) 

P  (^g) 

t  (sec) 

5  (km) 

ODE 

500.0 

85.0 

10.0 

777.41 

93.388 

6.705898 

10.0150 

RRA 

500.0 

85.0 

10.0 

777.48 

93.389 

6.705893 

10.0150 

ODE 

500.0 

95.0 

10.0 

222.59 

86.612 

6.705898 

10.0150 

RRA 

500.0 

95.0 

10.0 

222.52 

86.611 

6.705893 

10.0150 

ODE 

500.0 

135.0 

10.0 

545.01 

134.997 

9.448814 

14.1104 

RRA 

500.0 

135.0 

10.0 

545.00 

134.997 

9.448829 

14.1104 

For  SSP  4  (parabolic  SSP)  with  vertex  c(500)  =  1490  m/s,  speed  of  sound  at  the 
ocean  surface  c(0)  =  1500  m/s,  and  speed  of  sound  at  the  ocean  bottom  c(1000)  =  1500 
m/s,  i.e.,  an  ocean  depth  of  1000  m;  with  c(y)  =  c(500)  +  a(y  -  500)2  where 
a  =  0.00004,  at  a  horizontal  range  r—  10  km,  the  differences  between  the  two  algo- 
rithms are  as  follows: 

•  between  0.01  and  0.07  m  in  depth  y, 

•  0.001  deg  in  angle  of  arrival  /?  and, 

•  between  5  and  15  \x  sec  in  travel  time  t. 

Table  7,  when  compared  with  Table  4,  both  pertaining  to  an  arc  length  step  size  of 
2.0  m,  demonstrates  that  a  more  realistic  gradient  of  0.02  1/sec  enables  the  RRA  algo- 
rithm to  handle  propagation  in  a  SOFAR  channel  to  a  range  of  10  km  with  a  much 
improved  degree  of  accuracy.  Of  course,  this  could  be  improved  still  further  by  reducing 
the  arc  length  step  size. 
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B.     RRA  THREE-DIMENSIONAL  CAPABILITIES 

Figure  5  demonstrates  the  ability  of  the  RRA  algorithm  to  handle  launch  angles 
other  than  4>0  =  90°  .   This  figure  was  produced  using  the  parabolic  SSP  4  with 

c(y)  =  c{y0)  +  a(y  -  v0)2  ,  0  <  y  <  200  m,  (3.1) 

where  r0  =  100  m,  c{v0)  =  1490.0  m/s,  a  =  0.001  -^~  ,  and  </>„  =  85°  .  Note  that  the 
ray  stays  in  the  plane  created  by  this  angle,  as  illustrated  by  the  projection  of  the  ray 
onto  the  X  Z  plane  as  a  straight  line. 

The  RRA  algorithm  also  has  the  capability  to  divide  the  ocean  into  separate  zones 
as  a  function  of  z.  Each  zone  can  have  its  own  SSP  as  a  function  of  both  cross  range  jc 
and  depth  y.  For  computer  simulation  purposes,  the  speed  of  sound  was  defined  as  fol- 
lows: 

c{x,  y)  =  c(y0)  +  a(y  - y0)2  +  gx  x  ,  0  <  y  <  200  m  ,  (3.2) 

where^0  =   100  m,  c(y0)  =  1490.0  m/s,  a  =  0.001    -^j-  and 

f        0.0    1/s,    i  <  5  km 
x      [-0.05    1/s,  5  <  z  <  10  km 

The  particular  gradient  in  the  x  direction  was  selected  solely  for  computer  simulation 
purposes  to  illustrate  the  capability  of  the  RRA  algorithm  to  accommodate  sound  speed 
as  a  function  of  both  jc  and  y.  Figure  6  illustrates  the  ray  trace  provided  by  the  RRA 
algorithm  for  the  SSP  given  by  Equations  3.2  and  3.3,  showing  the  ray  staying  in  the 
Y  Z  plane  until  reaching  5  km,  where  the  gradient  in  the  x  direction  comes  into  effect, 
resulting  in  out-of-plane  propagation. 
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IV.     REPRESENTATION  OF  SOUND-SPEED  DATA 

A.     COMPARATIVE  ANALYSIS  -  SPATIAL  FOURIER  SERIES  VERSUS  AKIMA 
CUBIC  SPLINE  REPRESENTATIONS  OF  SOUND-SPEED  DATA 

Five  difFerent  exact  mathematical  functions  were  used  to  generate  1 1  sound-speed 
data  points  evenly  spaced  in  depth  in  a  200  m  deep  ocean.  Akima  cubic  spline  and 
spatial  Fourier  series  representations  of  the  SSPs  were  produced  from  these  data  points. 
Percent  error  calculations  of  their  sound-speed  fits  as  well  as  their  first  and  second-order 
derivative  fits  were  performed  with  respect  to  the  known  exact  values.  If  the  error  be- 
tween the  exact  and  modeled  values  was  less  than  10~6  ,  percent  error  was  set  to  zero. 

The  exact  functions  examined  were  as  follows: 


Constant  speed  of  sound  (  c  =   1500  m's  ) 

Linear  SSP  with  a  positive  gradient  (  c(y)  =   1500.0  +  0.017^  m/s  ) 

Linear  SSP  with  a  negative  gradient  (  c(>)  =  1500.0  —  0.01 7 y  m/s  ) 

rev 
Half-period  sine  wave  SSP  (  c(y)  =  1500.0  -  25.0  sin(  -~r  )  m/s  ) 

Parabolic  SSP  (  c(y)  =  1450.0  +  0.005(V  -  100.0)'  m/s  ) 


Tables  8  and  9  with  Figures  7,  8  and  9  illustrate  the  difficulties  that  the  Fourier 
representation  encounters  in  dealing  with  data  derived  from  an  exact  constant  speed  of 
sound.  No  noticeable  difference  in  the  sound-speed  ,/it  itself  is  observed  for  either  the 
Akima  cubic  spline  or  the  Fourier  series  representation,  but  the  latter  generates  some 
small  errors  in  both  the  first  and  second-order  derivatives. 

Tables  10  and  11  along  with  Figures  10,  11  and  12  compare  the  Akima  cubic  spline 
and  Fourier  series  representations  of  sound-speed  data  generated  from  a  linear  SSP  with 
a  positive  gradient.  In  Figure  10  and  Table  1 1  it  can  be  seen  that  the  P'ourier  represen- 
tation does  not  precisely  match  the  sound-speed  function,  although  Figure  10  and  Table 
10  show  that  the  Akima  cubic  spline  does.  Problems  similar  to  those  noted  for  a  con- 
stant speed  of  sound  can  be  observed  by  examining  Figures  11  and  12  and  Table  11. 
Discrepancies  in  the  derivatives  are  now  much  larger. 

Tables  12  and  13  along  with  Figures  13,  14  and  15  show  analogous  results  for  a 
linear  SSP  with  a  negative  gradient.  Once  again,  neither  the  sound-speed  fit  nor  the 
derivatives  are  precisely  modeled  by  the  Fourier  series  representation. 
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Table  8.     AKIMA    CUBIC    SPLINE     RESULTS    -    CON- 
STANT SPEED  OF  SOUND 
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Table  9.     FOURIER    SERIES 
SPEED  OF  SOUND 


RESULTS    -    CONSTANT 


o 

IX 

a: 

UJ 

.v 


o  - 

Q 

u 


ooooooooooooooooooooo 
ooooooooooooooooooooo 
ooooooooooooooooooooo 

ooooooooooooooooooooo 


-a-*>a--i-vj'-**>s-<j-v»-<i->r>j-'*^>j-vS->*v*-v*->3' 


I 


I 


I  I 


I  I 


I  I  I  I  I 


I  I 


III    ill    ni    in    11)    Hi    Ml    til    |]|    111    in    in    I    |||    111    hi    m    iii    \,f    t|i 

HOOH(MNOHNlflMNWe>NO>HNNMIfl 
U)UlONC0iT^ONN^N<]'Nlf|in>SHHNi)' 
UIO'SHKtMNHOIIIIOit^OI/lirHO^rilO 
NvOH^NNHNNrtrtHHHNHMHIflHxO 

ooooooooooooooooooooo 
III      II      II      II      III 


o 

tt 


o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

+  ♦  +  +  +  ♦♦ 

lit  ni  til  in  ill  til  ill 

o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

o  o  o  o  o  o  o 


o  o  o  o  o 

o  o  o  o  o 

+  ♦  +  ♦  ♦ 

tt  I  lit  It  I  lit  III 

o  o  o  o  o 

o  o  o  o  o 

o  o  o  o  o 

o  o  o  o  o 


ooooooooo 
ooooooooo 
+  ♦  +  +  ♦♦♦♦♦ 

HJHJHJHJHJHJHJHJUJ 

ooooooooo 
ooooooooo 
ooooooooo 
ooooooooo 


ooooooooooooooooooooo 


ooooooooooooooooooooo 
ooooooooooooooooooooo 
ooooooooooooooooooooo 

ooooooooooooooooooooo 


M'}'MKimiflMMlflMi},'tMMtfMKlMMMN> 


1- 
l-l 

u. 

„ 

1 

UJ 

1 

LU 

1     1 

UJ   UJ 

1 

UJ 

i 

UJ 

i 

UJ 

i 

UJ 

1 

UJ 

1 

UJ 

1     1     1     1     1     1     1     1     1     1     1 

11  1     1  l  1     Hi     HI     III     III     III     l||     I i I     LU     111 

w 

Ul 

vO 

un 

O*    r-t 

nD 

(^ 

CM 

to 

CM 

-o 

(O^o^m^ioNvoff'ifl 

h- 

\ 

nJ" 

cr 

-s-  un 

un 

CO 

CM 

o 

<r 

a 

a-cocOvj-uncor-^^-r^CMf^- 

o 

•H 

m 

un 

-*  -3- 

* 

CM 

H\ 

f-l 

s± 

■& 

MN^MfllflNlflOMH 

a 
u 

— 

vO 

r~- 

^   CM 

r-l 

.H 

i-H 

l-l 

r-i 

r-l 

if'^HrlCOHHNNvSIO 

o 

o 

O   O 
1       1 

O 

o 

o 
1 

o 
1 

o 

o 

ooooooooooo 
II               II               III 

_ , 

o 

o 

O  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

►- 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

( > 

^ 

+ 

+ 

+■    ♦ 

■*• 

+ 

♦ 

*■ 

♦ 

+ 

+    ♦♦♦♦♦♦    +    ♦    +    ♦ 

< 

o 

Ul 

UJ 

UJ   UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

til    III   til   III    III   III   III   lit   III   III    til 

X 

s 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

UJ 

iH 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

w 

w 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

t- 

l- 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

t-t 

M 

o 

r* 

a 

o 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

II 

UJ 
M 

rr 

Of 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

i- 

UJ 

o 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

r  V) 

CE 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

H 

CE 

O 

a 

UJ 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

a. 

UJ 
M 

x- 

Q 

< 

a: 

> 

«* 

o 

*" » 

^ 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

£> 

o 

u- 

1- 

Ul 

</> 

l-l 

s 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

—i 

_i 

u. 

r 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

u_ 

«* 

< 

w 

un 

un 

un  un 

un 

un 

un 

Ul 

un 

un 

ununununununununununui 

O 

z 

M 

M 

1- 

u 

o 

CD 

< 

■^ 

Ul 

H 

a 

t- 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

UJ 

tt 

01 

o 

^ 

DL 

o 

$ 

ul 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

«/) 

•• 

£ 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

u. 

UJ 

UJ 

un 

un 

un  un 

un 

un 

un 

un 

un 

un 

ununununununu-iunununun 

t 

o 

3 

u 

^ 

tt 

M 

1- 

UJ 

t: 

X 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

i 

(fl 

X 

t- 

— • 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

o 

a. 

T. 

o 

2 

UJ 

UJ 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

ooooooooooo 

u 

»- 

a 

!-( 

CM  M 

* 

un 

.0 

c~ 

CO 

<r 

OHRIM^lfl^NlOO^O 

26 


\<r>\ 

,uj| 

uj'EI 

d 

BI|**I 

a* 

B 

i.i 

•HlO 

OOll^l 

yn 

n 

21 

U) 

o,_i 

E 

i 

.1 

a: 

y 

<C  — 

1- 

«ia:l 

s 

i  > 

vil; 

Ul 

<r,cr> 

• 

■ 

a 


en 
_J 
o 

Q_ 


a= 
o 

LJ 


o 

CO 
LJ 


o 

LJ 




j                 i                 i 

!                 :                 1 
ill: 

1 1 1 1 1 

i                       i                       ;                       i 

!               j               !               j 
j               |               j.              | 

|                    i                    i                    i 

i                     i                      i                      i 
!                      1                      !                      i 
j                     j                     {                     j 

I                    !                    i                    ! 

III! 

'             *             7             T             T             1 
■                 I                 J                 i                 i 

7             7             7             7 

1                j                j                i 
:                 j                 :                 ! 



;j|; 

j                j                j                j                j                ;                ■                |                j 

j                j                |                j                i                i                j                j                j 

8 

i                 i                 i                 i 
9              I              8               i 

H               9               B               8               | 

-3 


§£ 


5 

S.   « 

z 

tar     ■— 
i  e 

poo 

in  o 

U.       O       ft 

o 

D£  en 

a      u  O 

Q     <r>  rv 


8     § 


in     a»      x 


iu)Hi<aa 


Figure  7.      Sound  Speed  Fit  -  Constant  Speed  of  Sound. 
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Figure  8.      First-Order  Derivative  Fit  -  Constant  Speed  of  Sound. 
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Figure  9.      Second-Order  Derivative  Fit  -  Constant  Speed  of  Sound. 
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Figure   10.      Sound  Speed  Fit  -  Linear  SSP  with  a  Positive  Gradient. 
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Figure  11.      First-Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient. 


33 


!   1 

'en 

<u 

ia: 

t»J,LJ 

d 

=c|«o| 

9$ 

SlS 

Q 

iae 

a  o 

UiD 

—  o: 

§H 

<n 

CJIJ 

CCl— I 

JT,I-. 

g 

x¥\ 

<E  (Ol 

o 

O 

o 

e— 
a: 

_] 
o 

Q_ 
LJ 


en 

t — I 

en 


a: 
o 


en 


(uiKtcoa 


Figure  12.      Second  Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient. 
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Table   13. 


FOURIER  SERIES   RESULTS  -   LINEAR 
WITH  A  NEGATIVE  GRADIENT 
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Figure  13.      Sound  Speed  Fit  -  Linear  SSP  with  a  Negative  Gradient. 
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Figure  14.      First-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient. 
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Figure  15.      Second-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient. 
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The  above  three  exact  sound-speed  functions  were  used  to  compare  the  Fourier  se- 
ries representation  with  the  Akima  cubic  spline  since  the  Fourier  fit  can  be  expected  to 
be  at  its  worst  since  we  are  trying  to  form  linear  functions  from  trigonometric  ones.  The 
next  two  functions  examined,  however,  are  non-linear;  the  half-period  sine  wave  and  the 
parabola. 

Tables  14  and  15  along  with  Figures  16,  17  and  18  compare  the  Akima  cubic  spline 
and  Fourier  series  representations  of  sound-speed  data  generated  from  a  half-period  sine 
wave.  Both  techniques  encounter  some  difficulty  here,  with  the  Fourier  representation 
generally  behaving  better  near  the  midpoint  of  the  data  than  the  Akima  cubic  spline  fit 
for  sound  speed  and  both  derivatives.  However,  near  the  surface  and  bottom,  the  errors 
in  the  Fourier  representation  become  much  larger  than  those  of  the  Akima  cubic  spline 
fit. 

Tables  16  and  17  along  with  Figures  19,  20  and  21  compare  the  Akima  cubic  spline 
and  Fourier  series  representations  of  sound-speed  data  generated  from  a  parabola. 
Surprisingly,  the  Akima  cubic  spline  representation  matches  sound  speed  and  derivatives 
very  well.  The  Fourier  representation  behaves  in  a  manner  similar  to  its  handling  of  the 
half-  period  sine  wave.  It  produces  small  errors  near  the  midpoint  of  the  data  and  large 
discrepancies  at  either  end  for  sound  speed  as  well  as  in  the  first  and  second-order  de- 
rivatives. 

The  disappointing  performance  of  the  Fourier  series  representation  of  sound  speed 
and  of  the  first  and  second-order  derivatives  of  sound  speed  led  to  various  attempts  to 
salvage  it.  First,  Lanczos  smoothing  as  discussed  in  Chapter  II,  Section  C,  was  imple- 
mented. This  ameliorated  the  oscillations  of  the  Fourier  representation  as  shown  in 
Tables  18  and  19  and  Figures  22  through  27  for  the  half-period  sine  wave  and  parabolic 
SSPs.  It  can  be  seen  that  Lanczos  smoothing  results  in  a  smoother  fit  for  sound  speed 
and  for  the  derivatives  to  a  point  closer  to  the  end  points  of  the  data  but,  nevertheless, 
leaves  significant  discrepancies  near  the  end  points. 
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Table  14.     AKIMA       CUBIC       SPLINE       RESULTS 
HALF-PERIOD  SINE  WAVE  SSP 
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Table   15. 


FOURIER  SERIES  RESULTS  -  HALF-PERIOD 
SINE  VVAVESSP 
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Figure  16.      Sound  Speed  Fit  -  Half-Period  Sine  Wave  SSP. 
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Figure   17.      First-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP. 
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Figure  18.      Second-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP. 
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Table  16.     AKIMA       CUBIC       SPLINE       RESULTS 
PARABOLIC  SSP 
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Table   17. 


FOURIER   SERIES    RESULTS 
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Figure   19.      Sound  Speed  Fit  -  Parabolic  SSP. 
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Figure  20.      First-Order  Derivative  Fit  -  Parabolic  SSP. 
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Figure  21.      Second-Order  Derivative  Fit  -  Parabolic  SSP. 
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Table   19.     FOURIER  SERIES  RESULTS  WITH  LANCZOS 
. SMOOTHING  -  PARABOLIC  SSP      , 
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Figure  22.      Sound  Speed  Fit  uith  Lanczos  Smoothing  -     Half- Period  Sine  Wave 
SSP. 
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Figure  23.      First-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Hall- Period  Sine 
Wave  SSP. 
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Figure  24.      Second- Order  Derivative  Fit  with  Lanczos  Smoothing  -  Half-Period  Sine 
Wave  SSP. 
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Figure  25.      Sound  Speed  Fit  with  Lanczos  Smoothing  -  Parabolic  SSP. 
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Figure  26.      First-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Parabolic  SSP. 
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Figure  27.      Second-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Parabolic  SSP 
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The  difficulty  in  achieving  a  good  fit  at  or  near  the  end  points  of  the  data  using  a 
spatial  Fourier  series  representation  is  a  result  of  the  instantaneous  commencement  and 
cessation  of  sound-speed  data  at  the  ocean  surface  and  bottom.  Therefore,  an  attempt 
to  overcome  this  difficulty  was  made  by  adding  synthetic  data  points  above  the  ocean 
surface  and  below  the  ocean  bottom  so  as  to  enable  the  Fourier  series  to  "settle  down" 
before  the  real  data  begins. 

Synthetic  data  was  created  by  first  examining  the  three  data  points  nearest  the  ocean 
surface  and  bottom.  Similar  procedures  were  used  to  generate  synthetic  data  both  above 
the  surface  and  below  the  bottom.   Therefore,  only  the  surface  case  will  be  examined. 

Let 


c(y{)_-  c(0) 

y\ 


"*io=  ,.  (4-1) 


where  c(0)  =  c{y  =  0)  ,  and 


cjv2)  ~  <j>i) 

J'2    "  >'l 


"21-        A-V  <4"2> 


where  j>,  is  depth  in  meters  at  the  ith  data  point.    If  we  let 

Am  =  m2\  —  m\o  (4-3) 

and  if 

| Am |    =    |m21   -  m]0\    <£  (4.4) 

for  some  arbitrarily  small  e,  then  linear  extrapolation  is  used,  i.e., 

c(y-i)  =  c(0)  +  mXQy_t,  /-  l,2,...,SYNPTS  (4.5) 

where 

y-t=  -yi  (4-6) 

and  SYNPTS  is  the  number  of  synthetic  data  points  desired. 
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If  the  test  at  Equation  4.4  fails,  and  if 

I  "Ho  I   <    l«2i  I  .  (4.7) 

then  mirror  imaging  about  the  ocean  surface  is  appropriate,  i.e., 

clyj  =  ctyd,  i  =  l,2,..,SYNPTS.  (4.8) 

If  both  of  the  tests  at  Equations  4.4  and  4.7  fail,  then  neither  linear  extrapolation 
nor  mirror  imaging  is  appropriate.  Under  these  circumstances  it  is  appropriate  to  re- 
verse the  rate  of  change  of  the  slope  and  let 

c{yj  =  c(0)  +  (m10  -  iAm)y_i,  i  =   \,2,.~,SYNPTS  .  (4.9) 

This  methodology  of  extending  the  real  data  was  implemented  using  five  synthetic 
data  points  both  above  the  surface  and  below  the  bottom.  The  results  are  shown  in 
Tables  20  through  24  and  Figures  28  through  42. 

A  substantial  improvement  in  the  accuracy  of  the  Fourier  series  representation  can 
be  seen  with  respect  to  all  the  functions  examined.  Nevertheless,  some  difficulty  is  still 
evident  near  the  ends  of  the  original  data.  A  final  attempt  to  resolve  this  discrepancy 
was  made  by  using  the  Akima  cubic  spline  representation  of  sound  speed  as  input  to  the 
Fourier  algorithm,  and  sampling  the  SSP  generated  by  the  Akima  cubic  spline  fit  at  some 
arbitrarily  high  frequency.  This  was  attempted  using  a  sampling  rate  ten  times  greater 
than  available  from  the  original  data.  The  synthetic  data  as  described  above  was  not  fit 
using  the  Akima  cubic  spline  method.  Results  for  the  half-period  sine  wave  and 
parabolic  SSPs  are  presented  in  Tables  25  and  26  and  Figures  43  to  48. 

These  results  are  once  again  an  improvement  on  previous  work.  Nonetheless,  some 
inaccuracies  remain.  Unfortunately,  time  constraints  precluded  further,  more  sophisti- 
cated attempts  to  further  improve  the  accuracy  of  the  Fourier  representation. 
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Table  20.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  - 
CONSTANT  SPEED  OF  SOUND 
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Table  21. 


FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  -  LIN- 
EAR SSP  WITH  POSITIVE  GRADIENT 
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Table  22.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  -  LIN- 
EAR SSP  WITH  NEGATIVE  GRADIENT 
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Table  23. 


FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  - 
HALF-PERIOD  SINE  WAVE  SSP 
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Table  24.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  - 
PARABOLIC  SSP 
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Figure  28.      Sound-Speed  Fit  -  Constant  Speed  of  Sound  with  Synthetic  Data. 
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Figure  29.       First-Order  Derivative  Fit  -  Constant  Speed  of  Sound  with  Synthetic 
Data. 
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Figure  30.      Second-Order  Derivative  Fit  -  Constant  Speed  of  Sound  with  Synthetic 
Data. 
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Figure  31.      Sound-Speed  Fit  -  Linear  SSP  with  a  Positive  Gradient  and  Synthetic 
Data. 
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Figure  32.       First-Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient  and 
Synthetic  Data. 


70 


gi 


ja- 
il 


o 

a 

o 


o 


cc 
o 


LJ 


o 
l_ 

en 

LJ 


O 
LJ 


-t I-  — 


> 1- -™».. 


+ - I — i -» » + > t- 


t 


A 


r\ 


'  ...J..i/[ 


a 


M 


1/ 


a 


\ 


1 


V/l 


a 


t t 


1 


ft— t 


1 1 —4 


i M- 


l    I 


(wiKicoa 


-a       " 


S   ~ 


yes 

«     a     — 

_I  X 

l.     o     o 
o  o 

u      U      in 
>     m 
—     ^     tn 

£28 


Figure  33.      Second-Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient  anc 
Synthetic  Data. 
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Figure  34.      Sound-Speed  Fit  -  Linear  SSP  with  a  Negative  Gradient  and  Synthetic 
Data. 
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Figure  35.       First-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient  and 
Synthetic  Data. 
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Figure  36.      Second-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient  and 
Synthetic  Data. 
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Figure  37.      Sound-Speed  Fit  -  Half-Period  Sine  Wave  SSP  with  Synthetic  Data. 
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Figure  38.      First-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP  with  Synthetic 
Data. 
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Figure  39.      Second-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP  with  Syn- 
thetic Data. 
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Figure  40.      Sound-Speed  Fit  -  Parabolic  SSP  with  Synthetic  Data. 
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Figure  41.      First-Order  Derivative  Fit  -  Parabolic  SSP  with  Synthetic  Data. 
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Figure  42.      Second-Order  Derivative  Fit  -  Parabolic  SSP  with  Synthetic  Data. 
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FOURIER  SERIES  RESULTS  WITH  LANCZOS  SMOOTHING, 
SYNTHETIC  DATA,  AKIMA  CUBIC  SPLINE  INPUT  TO  FOURIER 
METHOD  -  PARABOLIC  SSP 
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Figure  43.      Sound-Speed  Fit  -  Half-Period  Sine  Wave  SSP,  Akima  Cubic  Spline 
Input  to  Fourier  Method. 
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Figure  44.      First-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP,  Akima  Cubic 
Spline  Input  to  Fourier  Method. 
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rigure  45.      Second-Order  Derivative  Fit  -  Half- Period  Sine  Wave  SSP,  Akima  Cu- 
bic Spline  Input  to  Fourier  Method. 
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Figure  46.      Sound-Speed  Fit  -  Parabolic  SSP,  Akima  Cubic  Spline  Input  to  Fourier 
Method. 
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igure  47.      First-Order  Derivative  Fit  -  Parabolic  SSP,  Akima  Cubic  Spline  Input 
to  Fourier  Method. 
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igure  48.      Second-Order  Derivative  Fit  -  Parabolic  SSP,  Akima  Cubic  Spline  Input 
to  Fourier  Metbod. 


B.     RAY  TRACE  RESULTS  USING  SPATIAL  FOURIER  SERIES 
REPRESENTATION  OF  SOUND-SPEED  DATA  WITH  THE  RRA  ALGORITHM 

1.     Sound  Speed  Including  Deterministic  Component  Only 

The  RRA  algorithm  using  the  simple  Fourier  series  sound-speed  representation 
including  Lanczos  smoothing  but  without  synthetic  data  or  Akima  cubic  spline  input  to 
the  Fourier  method  was  tested  using  the  same  source  locations,  SSPs,  and  rays  as  in 
Chapter  III,  Section  A.   That  is,  the  SSPs  were  once  again  defined  as  follows: 

•  SSP  1:  Constant  speed  of  sound  (c  =   1500a/z/.s) 

•  SSP  2:  Linear  SSP  with  a  positive  gradient  (  c(y)  =  c(0)  +  gy    m's,  where  the 
speed  of  sound  at  the  surface  c(0)  =  ciy  =  0)  =  1500  m/s  and  g  =    4-0.017    1/s  ) 

•  SSP  3:  Linear  SSP  with  a  negative  gradient  (  c(y)  =  c(0)  +  gy    m/s,  where  the 
speed  of  sound  at  the  surface  c(0)  =  c(y  =  0)  =  1500  m/s  and  g  =    —0.017    1/s  ) 

•  SSP  4:  Parabolic  SSP  with  vertex  at  y0  where  c(}'0)  =  1490  m/s  and  the  speed  of 
sound    at    both    the    ocean    surface    and    bottom    is    equal    to     1500    m/s 

(c(y)  =  c(y0)  +  a(y  -  100)2    m/s  with  a  =  0.001    -^j  for  an  ocean  depth  of  200 
m). 

The  sound  source  spatial  location  was  established  at  .v0  =  0  m  and  z0  =  0  m.  A 
source  depth  ofy0  =  100  m  was  established  for  an  ocean  depth  of  200  m,  unless  otherwise 
noted.  Launch  angles  /?„  were  45°,  85°  and  135°  for  SSPs  1,  2  and  3;  and  85°,  95°,  and 
135°  for  SSP  4,  which  again  was  designed  to  simulate  a  SOFAR  channel. 

The  ODE  algorithm  using  the  Akima  cubic  spline  fit  to  sound-speed  data  re- 
mained as  the  benchmark  for  comparison  purposes,  as  detailed  in  Chapter  III,  Section 
A.  For  ease  of  reference,  the  results  of  the  ODE  algorithm  are  reproduced  in  the  tables 
of  this  sub-section.  The  ODE's  range  step  size  was  left  at  2.0  m  and  the  arc  length  step 
size  in  the  RRA  algorithm  was  left  at  2.0  m  as  well,  except  where  otherwise  noted. 

The  results  of  the  various  computer  simulation  runs  are  presented  in  Tables  27 
to  33  below. 
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Table  27.     RAY  ACOUSTICS   RESULTS  FOR  SSP   1,  RAYS   1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

Vo  (m) 

fi>  (deg) 

r  (km) 

y  (m) 

P  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

100.00 

45.000 

9.428090 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.94 

45.000 

9.42S0S8 

14.1421 

ODE 

100.0 

85.0 

10.0 

174.89 

85.000 

6.692132 

10.0382 

RRA 

100.0 

S5.0 

10.0 

174.88 

S5.000 

6.692132 

10.0382 

ODE 

100.0 

135.0 

10.0 

100.00 

135.000 

9.428090 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.06 

135.000 

9.428088 

14.1421 

For  SSP  1  (constant  speed  of  sound),  Table  27  shows  that  at  a  horizontal  range 
of  r=  10  km,  the  differences  between  the  benchmark  ODE  algorithm  and  the  RRA  al- 
gorithm using  the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.01  and  0.06  m  in  depths  and 

•  2  fj.  sec  in  travel  time  t. 

Note  the  significant  degradation  of  the  RRA  results  relative  to  those  obtained 
using  the  Akima  cubic  spline  fit  as  in  Table  1  of  Chapter  III,  Section  A.  This  compar- 
ative degradation  of  consistency  between  the  ODE  and  RRA  algorithms  will  be  evident 
throughout  this  subsection. 
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Table  28.     RAY  ACOUSTICS  RESULTS  FOR  SSP  2,  RAYS   1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

y0  (m) 

Po  (cleg) 

r  (km) 

y(m) 

P  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.417413 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.97 

45.000 

9.417429 

14.1421 

ODE 

100.0 

85.0 

10.0 

162.5G 

85.487 

6.683492 

10.0374 

RRA 

100.0 

85.0 

10.0 

163.31 

85.482 

6.683550 

10.0374 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.417413 

14.1421 

RRA 

100.0 

85.0 

10.0 

100.03 

135.000 

9.417430 

14.1422 

For  SSP  2  (linear  SSP  with  a  positive  gradient),  Table  28  shows  that  at  a  hori- 
zontal range  r  =  10  km,  the  differences  between  the  benchmark  ODE  algorithm  and  the 
RRA  algorithm  using  the  Fourier  series  sound-speed  Fit  are  as  follows: 

•  between  0.02  and  0.19  m  in  depth  y, 

•  0.005  deg  in  angle  of  arrival  /?, 

•  between  16  and  58  ix  sec  in  travel  time  t,  and 

•  0.1  m  in  path  length  s. 
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Table  29.     RAY  ACOUSTICS  RESULTS  FOR  SSP  3,  RAYS   1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

>'o  (m) 

k  (deg) 

r  (km) 

y  (m) 

P  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.438784 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.96 

45.000 

9.438799 

14.1422 

ODE 

100.0 

85.0 

10.0 

168.89 

84.512 

6.699404 

10.0380 

RRA 

100.0 

85.0 

10.0 

169.76 

84.492 

6.699471 

10.0380 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.438784 

14.1421 

RRA 

100.0 

85.0 

10.0 

100.04 

135.000 

9.438798 

14.1421 

For  SSP  3  (linear  SSP  with  a  negative  gradient),  Table  29  shows  that  at  a  hori- 
zontal range  r=  10  km,  the  difference  between  the  benchmark  ODE  algorithm  and  the 
RRA  algorithm  using  the  Fourier  sound-speed  Fit  are  as  follows: 

•  between  0.03  and  0.13  m  in  final  depths, 

•  0.020  deg  in  angle  of  arrival  /?, 

•  between  14  and  67  /x  sec  in  travel  time  t,  and 

•  0.1  m  in  path  length  s. 

The  results  contained  in  the  above  tables,  while  poorer  than  those  detailed  in 
Chapter  III,  Section  A,  nevertheless  produce  travel  times  which  are  consistent  in  the 
millisecond  region.   However,  this  is  not  the  case  for  the  SSPs  examined  next. 
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Table  30.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS   1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

y0  (m) 

ft (dcg) 

r  (km) 

y  (m) 

P  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

37.06 

87.244 

6.710387 

10.0183 

RRA 

100.0 

45.0 

10.0 

42.86 

86.674 

6.709733 

10.0185 

ODE 

100.0 

85.0 

10.0 

162.94 

92.756 

6.710387 

10.0183 

RRA 

100.0 

85.0 

10.0 

157.14, 

93.326 

6.709733 

10.0185 

ODE 

100.0 

135.0 

10.0 

144.97 

134.922 

9.448852 

14.1104 

RRA 

100.0 

85.0 

10.0 

143.98 

134.930 

9.448515 

14.1112 

For  SSP  4  (parabolic  SSP),  Table  30  shows  that  at  a  horizontal  range  of  r  =  10 
km,  the  differences  between  the  benchmark  ODE  algorithm  and  the  RRA  algorithm 
using  the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.99  and  5.8  m  in  final  depth  y, 

•  between  0.008  and  0.570  deg  in  angle  of  arrival  /?, 

•  between  337  and  654  jx  sec  in  travel  time  r,  and 

•  0.2  m  in  path  length  5. 

As  in  Chapter  III,  Section  A,  SSP  4  shows  significantly  less  consistency  be- 
tween the  two  algorithms  than  the  previous  SSPs  examined.  Again,  two  approaches 
were  taken  to  attempt  to  reduce  this  discrepancy.  First,  the  arc  length  step  size  was  re- 
duced from  2.0  to  1.0  m  (see  Tables  12  and  13)  and  second,  the  steepness  of  the  gradient 
was  reduced  from  0.1  to  0.02    1/sec  (see  Table  14). 
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Table  31.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  FOURIER  SOUND-SPEED  FIT 


Method 

>'o(m) 

P0  (deg) 

r  (km) 

y  (m) 

fi  (deg) 

r  (sec) 

5  (km) 

ODE 

100.0 

45.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RRA 

100.0 

45.0 

5.0 

68.79 

85.435 

3.354591 

5.0090 

ODE 

100.0 

85.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

131.21 

94.565 

3.354591 

5.0090 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

85.0 

5.0 

77.82 

45.017 

4.724162 

7.0554 

For  SSP  4  (parabolic  sound  speed),  Table  31  shows  that  at  a  horizontal  range 
of  r  =  5  km,  using  an  arc  length  step  size  of  2.0  m,  the  differences  between  the 
benchmark  ODE  algorithm  and  the  RRA  algorithm  using  the  Fourier  series  sound-speed 
fit  are  as  follows: 

•  between  0.32  and  4.54  m  in  final  depth  y, 

•  between  0.002  and  0.162  deg  in  angle  of  arrival  /?, 

•  between  188  and  249  /j.  sec  in  travel  time  t,  and 

•  0.2  m  in  path  length  s. 

The  data  in  Table  31  above  was  generated  in  order  to  enable  a  direct  compar- 
ison with  the  arc  length  step  size  set  to  1.0  m,  as  shown  in  Table  32  below.  The  re- 
duction in  final  range  from  10  to  5  km  was  necessitated  by  the  structure  of  the  RRA 
computer  code,  which  required  large  quantities  of  computer  storage  space  to  be  allo- 
cated, particularly  as  arc  length  step  size  was  reduced. 
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Table  32.     RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  1.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

>'o  (m) 

ft (deg) 

r  (km) 

>'(m) 

fi  (deg) 

t  (sec) 

s  (km) 

ODE 

100.0 

45.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RRA 

100.0 

45.0 

5.0 

68.86 

85.442 

3.354604 

5.0090 

ODE 

100.0 

85.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

131.14 

94.558 

3.354604 

5.0090 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

85.0 

5.0 

77.74 

45.019 

4.724126 

7.0554 

For  SSP  4  (parabolic  SSP),  Table  32  shows  that  at  a  horizontal  range  of  r  =  5 
km,  and  with  the  RRA  arc  length  step  size  set  to  1.0  m,  the  differences  between  the 
benchmark  ODE  algorithm  and  the  RRA  algorithm  using  the  Fourier  series  sound-speed 
fit  are  as  follows: 

•  between  0.24  and  4.61  m  in  final  depth  y, 

•  0.155  deg  in  angle  of  arrival  /?, 

•  between  175  and  285  fi  sec  in  travel  time  t,  and 

•  0.2  m  in  path  length  s. 


95 


Table  33.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  BOTTOM  DEPTH  1000  M,  FOURIER 
SOUND  SPEED  FIT 


Method 

}'o  (m) 

Po  (deg) 

r(km) 

y(m) 

fi  (deg) 

t  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

777.41 

93.388 

6.705898 

10.0150 

RRA 

100.0 

45.0 

10.0 

768.46 

93.522 

6.705420 

10.0152 

ODE 

100.0 

S5.0 

10.0 

222.59 

86.612 

6.705898 

10.0150 

RRA 

100.0 

85.0 

10.0 

231.54 

86.478 

6.705420 

10.0152 

ODE 

100.0 

135.0 

10.0 

545.01 

134.997 

9.448829 

14.1104 

RRA 

100.0 

85.0 

10.0 

544.58 

134.997 

9.448180 

14.1107 

For  SSP  4  (parabolic  speed  of  sound),  with  vertex  c(500)  =  1490  m/s  and  speed 
of  sound  at  the  ocean  surface  c(0)  =  1500  m/s  and  speed  of  sound  at  the  ocean  bottom 
c(1000)  =  1500  m/s,  i.e.,  an  ocean  depth  of  1000  m;  at  a  horizontal  range  of  r  =  10  km, 
the  differences  between  the  benchmark  ODE  algorithm  and  the  RRA  algorithm  using 
the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.43  and  8.95  m  in  final  depths, 

•  0.134  deg  in  angle  of  arrival /?, 

•  between  478  and  649  fj.  sec  in  travel  time  r,  and 

•  between  0.2  and  0.3  m  in  path  length  s. 

Table  33  shows  no  significant  improvement  over  Table  30  for  SSP  4  with  the 
original  extreme  gradient.  The  errors  associated  with  the  Fourier  sound  speed  repre- 
sentation have  resulted  in  this  lack  of  improvement,  whereas  substantial  improvement 
was  observed  between  these  two  cases  for  the  Akima  cubic  spline  fit  with  the  RRA  al- 
gorithm (see  Tables  4  and  7  in  Chapter  III,  Section  A). 
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2.     Sound  Speed  Including  Deterministic  and  Random  Components 

In  order  to  investigate  the  effects  of  randomness  in  the  SSP,  a  model,  of  the 
power  spectral  density  of  the  random  component  of  the  speed  of  sound  cR(y)  had  to  be 
proposed.  The  most  straightforward  model  of  band-limited  white  noise  was  used,  as 
given  by 

(So,    IjJI   <  BW 
S^]  =  to,      \fy  I   >  BW  (4'l0) 

where  BW  is  the  bandwidth,  and  S0  is  the  constant  magnitude  of  the  power  spectrum  at 
all  spatial  frequencies^  . 

As  a  consequence  of  the  assumption  of  zero  mean  white  noise,  the  variance  of 
cR{y)  is  given  by 


and,  as  a  result, 


o2Cr=2S0BW  (4.11) 


2 


A  value  for  the  variance  of  cR(y)  in  shallow  water  of  2.5  x  10~7  was  obtained 
from  Brekhovskikh  and  Lysanov  [Ref.  8:  p.  208  ].  The  simulation  program  was  run  with 
randomness  incorporated  as  above  and  for  SSP  1  (constant  speed  of  sound).  When  the 
results  were  compared  with  those  from  the  RRA  algorithm  previously  run  using  the 
Fourier  series  representation  of  sound-speed  data,  differences  of  a  maximum  of  0.1  m  in 
final  depth  y,  0.002  deg  in  angle  of  arrival  /?,  and  2  \x  sec  in  travel  time  r  were  observed. 
The  RRA  algorithm  was  then  run  with  SSP  4  (parabolic  SSP  in  a  200  m  deep  ocean) 
incorporating  randomness  and  differences  of  a  maximum  of  0.22  m  in  final  depth  y  , 
0.018  deg  in  angle  of  arrival  /?,  and  16  n  sec  in  travel  time  x  were  observed. 

The  variance  was  then  increased  to  determine  at  what  magnitude  it  caused  the 
randomness  to  become  significant.  Further  simulation  runs  revealed  that  for  SSP  1, 
differences  in  travel  times  of  greater  than  10  /j.  sec  needed  a\R  to  be  at  least  of  the  order 
of  10-3  ,  four  orders  of  magnitude  greater  than  the  estimate  given  in  Brekhovskikh  and 
Lysanov  [Ref.  8:  p.  208  ]. 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  recursive  ray  acoustics  (RRA)  algorithm  was  shown  to  be  an  accurate  ray 
acoustics  algorithm  when  used  in  conjunction  with  an  Akima  cubic  spline  fit  to  sound- 
speed  data.  The  RRA  algorithm  can  be  made  as  accurate  as  desired  by  reducing  the  arc 
length  step  size.  Although  the  ODE  algorithm  computes  amplitudes  which  the  RRA 
algorithm  does  not,  the  RRA  algorithm  was  observed  to  be  significantly  faster  in  terms 
of  computer  execution  time. 

The  ability  of  the  RRA  algorithm  to  handle  sound  speed  as  a  function  of  all  three 
spatial  coordinates  enables  its  application  to  real  world  problems  wherein  a  horizontal 
gradient  of  sound  speed  exists,  either  in  the  cross-range  or  down-range  direction.  This 
could  conceivably  be  the  case  for  either  long  range  ray  propagation  or  for  propagation 
in  estuaries  or  across  current  boundaries  where  significant  horizontal  gradients  can  exist. 

A  spatial  Fourier  series  (SFS)  representation  of  sound-speed  data,  while  attractive 
for  reasons  of  its  spectral  content  and  the  ease  of  incorporation  of  randomness  at  spe- 
cific harmonics,  encountered  difficulty  with  regard  to  providing  accurate  fits  to  sound- 
speed  data.  Several  techniques  were  used  to  attempt  to  rectify  this,  including  the  use 
of  Lanczos  smoothing  ,  addition  of  synthetic  data  above  the  ocean  surface  and  below 
the  ocean  bottom,  and  sampling  of  the  Akima  cubic  spline  fit  to  provide  data  points  for 
input  to  the  SFS  method.  These  techniques  resulted  in  some  substantial  success  in  the 
reduction  of  the  inaccuracies  of  the  SFS  method.  However,  although  the  SFS  method 
eventually  provided  accurate  fits  to  the  sound-speed  data,  it  still  provided  inaccurate  first 
and  second-order  derivatives  when  compared  against  the  results  obtained  from  the 
Akima  cubic  spline  method. 

A  comparison  was  made  between  the  ODE  algorithm  using  an  Akima  cubic  spline 
sound-speed  representation  and  the  RRA  algorithm  using  a  SFS  sound-speed  represen- 
tation, with  the  SFS  representation  in  its  most  basic  form,  that  is,  with  none  of  the 
techniques  referred  to  in  the  above  paragraph  incorporated.  The  RRA/SFS  combina- 
tion was  found  to  be  wanting  in  accuracy. 

As  was  mentioned  previously,  the  RRA  algorithm  using  the  Akima  cubic  spline  fit 
to  sound-speed  data  was  proven  fast  and  accurate,  and  can  be  made  as  accurate  as  de- 
sired by  reducing  the  arc  length  step  size.    However,  investigation  is  warranted  to  de- 
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termine  if  the  RRA  algorithm  can  be  generalized  to  include  amplitude  calculations, 
especially  at  focal  points  and,  hence,  caustics. 

Throughout  the  simulations  undertaken  in  this  thesis,  a  zero-slope,  flat  ocean  sur- 
face and  bottom  was  assumed.  In  order  to  apply  the  RRA  algorithm  to  realistic  ocean 
boundaries,  the  ability  to  accommodate  arbitrarily  shaped  surfaces  should  be  included. 
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